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ABSTRACT 

We present a new method to standardize Type la supernova (SN la) luminosities to S 0.13 magnitudes using flux ratios from a single 
flux-calibrated spectrum per SN. Using Nearby Supernova Factory spectrophotomery of 58 SNe la, we performed an unbiased search 
for flux ratios that correlate with SN la luminosity. After developing the method and selecting the best ratios from a training sample, 
we verified the results on a separate validation sample and with data from the literature. We identified multiple flux ratios whose 
correlations with luminosity are stronger than those of light curve shape and color, previously identified spectral feature ratios, or 
equivalent width measurements. In particular, the flux ratio 7?642/443 = F(642 nm)/F(443 nm) has a correlation of 0.95 with SN la 
absolute magnitudes. Using this single ratio as a correction factor produces a Hubble diagram with a residual scatter standard deviation 
of 0.1 25 ±0.01 1 mag, compared with 0.161+0.015 mag when fit with the SALT2 light curve shape and color parameters x\ and c. The 
ratio 7? 64 2/443 is an effective correction factor for both extrinsic dust reddening and instrinsic variations such as those of SN 1991T-like 
and SN 1999aa-like SNe. When combined with broad-band color measurements, spectral flux ratios can standardize SN la magnitudes 
to ~ 0. 12 mag. These are the first spectral metrics that give robust improvements over the standard normalization methods based upon 
light curve shape and color, and they provide among the lowest scatter Hubble diagrams ever published. 
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1. Introduction to correct observed magnitudes (e.g.. lRiess et all Fl998l) . While 

■ these methods showed promise for what could be acheived with 

\ The ability to standardize the absolute luminosity of Type la su- spectra, none of them is robustly competitive with light curve 

O ■ pernovae (SNe la) makes them powerful cosmological probes shape and color corrections alone. 

•• ■ through their use as luminosity distance indicators. Uncorrected In this work, we take a new unbiased approach to the prob- 
. ^ \ observations of SN la absolute magnitudes have an RMS scatter lem of finding optimal spectral flux ratios for standardizing 
K> , of -40%; corrections for light curve shape and color normalize SN la peak brightness. Flux ratios best correlated with absolute 
$_| ■ them to 15% to 20%, i.e., luminosity distances of SNe la can magnitudes are detected in a training set and cross-validated on 
\ be measure d to 7-10% accuracy using current standard i zation a separate subsample. Hubble diagram fits are compared using 
techniques dPhillipsL Il993t iPrieto et all l2006t iGuy et all l2007t standard lightcurve color and shape corrections, using only flux 
iJha et all |2007j). Improved standardization would reduce both ratios as standardization parameters, or using flux ratios corn- 
statistical and systematic errors, since the remaining dispersion bined with light curve color measurements. We find several new 
represents an upper limit to possible bias from uncorrected phys- robust spectral indicators that outperform light curve shape and 
ical effects that could arise from variations in host galaxy dust color corrections, 
properties, progenitor properties such as mass, metallicity, and 
age, and explosion physics such as geometry. 

Spectral indicator s such as feature rati o s {e.g. , iNugent et all 

1995; Bongar cfet all l2006t iFolev et all 120081) and pseudo- This analysis uses spectrophotometry of 58 of the supernovae 

equiv alent widths (e.g.. lHachinger et all 120061 iBronder et ail obtained by the Nearby Supernova Factory (SNfactory) collabo- 

2008) have been proposed to enhance peak SN la normaliza- ration using its Supe rNova Integral Field Spectrograph (SNIFS, 

tion from light curve shape and color, or even as alternative lAldering et all 120021) on the University of Hawaii 2.2-meter 

approaches altogether which use a single night's observations telescope on Mauna Kea. This subset consists of SNfactory 



2. Data and analysis 



SNe la for which final host-galaxy follow-up has been com- 
pleted and which pass minimal cuts on the quality of the light 
curve fit. Targets were additionally required to have a spectrum 
within +2.5 restframe days of B-band maximum light; if multi- 
ple spectra were available, the one nearest maximum light was 
used. Hubble residuals, colors, and light curve shape parame- 
ters were not a factor in selecting this subsample. These SNe are 
evenly distributed in the redshift range 0.02 < z < 0.09, which 
allows relative distance measurements with minimal uncertain- 
ties from peculiar velocities. The methods were developed on a 
training subset of 28 SNe before cross checking the results on 
a separate validation sample of 30 SNe to ensure that the final 
results are not simply overrating statistical fluctuations of the 
dataset. 

The subsets have consistent distributions of light curve width 
and color, redshifts, and Hubble residuals; these are consistent 
with other published SNe la used for cosmology measurements 
and training lightcurve fitters. The host galaxies include spirals, 
ellipticals, and irregulars, and the host properties span an order 
of magnitude in progenitor age, metallicity, and inferred amount 
of dust extinction. This is comparable to the estimated evolu- 
tion of the universal mean of these quan tities to z ~ 1 (e.g., 
iRiess & Livioll2006tlKowalski et aU l2008) thus the results from 
this dataset are also likely to work well for high-z SNe la. 

The spectra were dereddened to correct for Milky Way dust 
dCardelli et all 119891: ISchlegel et all 1 19981) . deredshifted, and 
rebinned in cAA/A ~ 2 000 km/s bins from 350 to 850 nm (23 to 
56 A bins) to equally sample spectral features in the physically 
relevant velocity space while still over-sampling SN la spectral 
features. The exact bin sizes used were not rigorously optimized 
for this analysis, though cross checks were performed to ensure 
that the results are not sensitive to the exact binning chosen. 

Photometry was synthesized from the multi-epoch flux- 
calibrated spectra in box filters corresponding to app roximately 
B, V, and R. Light curves were fit using SALT2 dGuv et all 
120071) to estimate their restframe B-band peak magnitudes mj 
and calculate their uncorrected B-band Hubble diagram resid- 
uals (or equivalently, their uncorrected absolute magnitudes). 
Alternatively, w# was synthesized directly from the spectra, in- 
dependent of an SN light curve model. 

The correlation of absolute magnitudes with each possible 
binned flux ratio was calculated for the training sample as shown 
in Fig. [T] The upper-left triangle of the color density plot shows 
the absolute Pearson correlation p of the flux ratios Hy/x - F y /F x 
with the SN la absolute magnitudes; the bottom-right triangle 
shows the color corrected correlations as described below. Using 
\ogCR) produced similar though slightly worse results. 

Spearman rank order correlations were also calculated in or- 
der to be less sensitive to outliers. Although the Spearman corre- 
lation is a more robust statistic, the Pearson correlation is a better 
predictor of the effects on a Hubble diagram ;f 2 fit since outliers 
affect the Pearson correlation and x 2 fits in similar ways. We 
selected the best ratios using the simple average of the Pearson 
and Spearman correlations; by considering both types of corre- 
lations, we select relatively strong correlations while avoiding 
those which are primarily due to outliers. We also checked the 
robustness of selected peaks by comparing the results from other 
near-max spectra in the training sample which were not the spec- 
tra closest to maximum light. 

The statistical significance of prospective correlations in the 
training set was assessed using Monte Carlo permutation testing 
before the validation set was unblinded. For each random trial, 
Hubble diagram residuals were randomized and the ratio with 
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Fig. 1: Absolute Pearson correlations of flux ratios H y/x = F y /F x with 
SN la absolute magnitudes. The upper-left triangle shows the correla- 
tions with uncorrected absolute magnitudes. For illustration, the lower- 
right triangle shows the correlations of color-corrected ratios K c with 
color-corrected magnitudes. A median SN la spectrum is included on 
each axis for reference along with the wavelengths used in the flux ra- 
tios of TableQ] 
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Fig. 2: Hubble diagram residuals vs. 7?642/443 for the training set (filled 
circles), validation set (open circles), and literature data (x symbols). 



the highest correlation was retained. We find that V{\p\ > 0.9) < 
10~ 5 , compared with multiple peaks above 0.9 in Fig.Q] 

After selecting the five best flux ratios from the training sam- 
ple, we unblinded the results for those ratios for the validation 
sample. Figure [2] shows the results for the highest ranked ra- 
tio, ^642/443- It plots the uncorrected Hubble diagram residuals 
AfiB vs. ??642/443 for the training set of supernovae (filled cir- 
cles) and validation set (open circles). The correlation is very 
strong (0.95), consistent between the two samples, and is not 
dominated by just a few supernovae. For comparison, the corre- 
lation of Hubble residuals with the combination of the SALT2 
light curve shape and color parameters x\ and c is 0.92. 

We fit a Hubble diagram of the 58 SNe using flux ratios % 
with a distance modulus of fig - (nig - M') + yR, where mg 
is the restframe B-band magnitude from either SALT2 (with- 
out any color or x\ corrections) or directly synthesized from 
the spectrophotometry, M' absorbs offsets from the % intercept 
and the mean SN la absolute magnitude, and y is solved in 
the fit along with M' to minimize the x 2 of the residuals. The 
result for 7^642/443 when using SALT2 is shown in Fig. [3] 
(top). The standard deviation of the residual scatter is <x(A/7g) = 



Table 1 : Correlations of flux ratios R with SN la absolute magnitudes and standard deviations of Hubble diagram fits. 
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Notes: 7 is a fit parameter in the distance modulus /j b = (m B - M') + yH, and cr core = 1.4826 x median(|A// s - median(A/j B )|). 
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Fig. 3: Hubble diagram residuals when fit with the single parameter 
7?642/443 (top), with color and dereddened ratio 7?g 41/443 (middle), and 
with the SALT2 parameters Xi and c (bottom). Gray bands highlight 
±lcr scatter. Intrinsically bright SN 1991T-like and SN 1999aa-like SNe 
are highlighted in black; highly dust-extincted SNe are highlighted with 
red squares. 



0.128 + 0.012 mag when using only 7?642/443 to standardize the 
peak magnitude mg, compared with 0.161 +0.015 mag when us- 
ing fig = (ntB - M) + ax\ - /3c instead (Fig. [3] bottom). In this 
case, a and /3 are allowed to vary in the fit. 

Synthesizing nig directly from the flux-calibrated near-max 
spectra provides a peak magnitude estimate independent of a 
light curve model. In this case ct{A^.b) - 0.125+0.01 1 mag when 
measuring mg and 7?642/443 from the same spectra. Note that this 
standardization is achieved using only a single flux calibrated 
spectrum per SN; the light curve fit is used only to determine the 
phase. This demonstrates the standardization power of a single 
spectrum and that these flux ratios are not merely correcting for 
some anomaly of the SALT2 light curve fit. 

In Fig.Ql intrinsically bright SN 1991T-like and SN 1999aa- 
like SNe are highlighted in black; highly extincted SNe are high- 
lighted with red squares. ^42/443 provides an effective correction 
for both of these extremes for which x\ and c are less effective. 

We also repeated this analysis fitting a combination of color 
and flux ratios. In this case, we dereddened the spectra us- 
ing the SALT2 color model so that the resulting spectra have 
E(B-V) = 0; this removes color from both extrinsic dust redden- 
ing as well as possible intrinsic color without trying to distinuish 
between the two effects. The remaining spectral variations, how- 
ever, should be predominantly due to small-scale intrinsic ef- 



fects and not extrinsic dust reddening. We denote flux ratios 
from these decolored spectra as K . Their correlation with color- 
corrected Hubble diagram residuals is shown in the lower-right 
triangle of Fig. [U the final selection was based upon the com- 
bined fit of c and 7? f , not just the correlations in Fig. [T] Several 
combinations of c and < R C have correlations with absolute magni- 
tude of ~ 0.95 and result in ct(A//b) ~ 0.12 mag. Figure [3] (mid- 
dle) shows the results for ^6 42/ 4 43 with cr(Afi B ) = 0.1 19 ± 0.01 1 
mag. Adding the light curve shape parameter x\ to the fits does 
not improve the results for the single-?? fits or the c, 7? c fits. 

The results are summarized in TablefT] Table|2]in the online 
version lists the supernovae used with their redshifts, x\, c, Afig, 
and f$ c) values. For all selected ratios, a single flux ratio results 
in smaller or competitive Hubble diagram residual scatter com- 
pared to what is acheived by standard corrections for light curve 
shape and color (x\ and c); the combination of a color-corrected 
ratio 9? c and c consistently outperforms x\ and c. The dispersion 
of the core is also improved; flux ratios are not merely correcting 
a few outliers. The largest discrepancy between the training and 
validation results is for ^41/417, but even in this case, the vali- 
dation set has a residual scatter of only 0.185 ± 0.025 mag. The 
fits are error weighted but the standard deviation is not, which 
allows some cases where the combined fit has a scatter either 
less than or greater than the scatter of both of the individual fits. 

3. Discussion 

The 130 to 330 nm wavelength baseline of the flux ratios means 
that they are effectively color measurements using very narrow 
band filters, thus they can correct for extrinsic dust extinction in a 
manner similar to standard E(B - V) corrections. In fact, 7?642/443 
has a correlation of 0.83 with E(B - V) and four of the five SNe 
with the largest ^642/443 values have host galaxy ISM features in 
their observed spectra, indicating extrinsic dust extinction. 

However, 7?642/443 also appears sensitive to intrinsic SNe 
variations that are correlated with absolute magnitude. There are 
other dim SNe in this sample with large ^42/443 and no spectral 
evidence for host extinction, and ^642/443 is able to standard- 
ize absolute magnitudes better than color alone: using c only, 
the Hubble residuals have a standard deviation of 0.215 ± 0.020 
mag. Additionally, the characteristic size of the correlation peaks 
in Fig. Q]are typical of the widths of SN features, indicating a 
sensitivity to instrinsic variations of the supernovae. 

PHOENIX calcu lations based upon the W7 explosion model 
dBongard et all l2008h show a strong correlation of 7?642/443 with 



the absolute fi-band magnitude, driven by large intrinsic vari- 
ations of Fe II and Fe III line blends around 443 nm. The re- 
gion redward of 640 nm has relatively little intrinsic variation 
in these models and serves as a normalization anchor for the 
ratio. The < 7?642/44 3 correlation also appears in a new survey 
of 2D SEDONA dKasen et all 120061) models spanning nickel 
mass 0.4-1.0 M Q with 30 different viewing angles per model; 
these generally reproduce the known Phillips brighter-broader 
correlation. The models exhibit a similar correlation between 
■^642/443 and absolute magnitude as we found, but the correla- 
tion is even stronger than the theoretically modeled Phillips re- 
lation (D. Kasen, private communication). More modeling work 
is needed to understand the underlying physics of our results. 

It may seem surprising that a single parameter could correct 
for the completely unrelated effects of SN la intrinsic variability 
and extrinsic dust reddening. But note that the slopes for a dust 
correction and any intrinsic correction will vary with different 
wavelength pairs selected for the flux ratio. Some wavelength 
pairs will align the slopes more than others, thus reinforcing the 
individual corrections and allowing for two physically distinct 
phenomena to be corrected with a single parameter. The final 
best choice will also depend upon the strength of the intrinsic 
correlations for that ratio and how accurately it can be measured. 

We checked our resu lts for ^47/443 u s ing sp ectra within 2.5 
days of maximum from Mathes on et alj J2008) after rejecting 
3 exhibiting strong H-a indicative of host galaxy contamina- 
tion. Uncertainties on 7?642/443 were assigned based on the inter- 
nal dispersion for the 12 SNe la having multipl e suitable spec- 
tra; ot herwise we used uncertainties derived by Mat heson et al.l 
(2008) from their comparison with photometr y. Corresponding 
.fi-b and peak magnitude s for 16 SNe la from iJha et al. I (12006b 
and lHicken et al.l d2009l) were put on a consistent scale and stan- 
dardized using 7?642/443, including photometric and redshift er- 
rors due to peculiar velocities. With the exception of SN 1999cl, 
we find good agreement with our results with a consistent slope 
y and strong correlation p = 0.93 (see Fig. |2). The scatter is 
cr(AfiB) = 0.25 mag (0.22 for an error weighted <x), but the fit 
has xl - 0.9, indicating that the dispersion is consistent with 
the measurement errors for these data and is not a contradiction 
of our low dispersion result. The environment o f SN 1999cl is 
know n to be unusual, having time-vary ing Na I D (Blondin et all 
2009) and "highly nonstandard dust" (Krisciuna s et all l 2006l) . It 
is also a large outlier (1.5 mag) when corrected by x\ andc. This 
exception is consistent with the inability of any single parameter 
to simultaneously correct intrinsic SN variations and more than 
one type of dust behavior. 

Application of ^642/443 requires flux calibrated spectra near 
maximum with accurate host galaxy subtraction; spectral slope 
errors should be kept below a few percent. To keep the noise 
contribution to the Hubble residuals under 0. 1 mag, S /N <: 25 
in 2 000 km/s bins is sufficient. For distant SNe la, where S /N 
may be a limiting constraint, the spectra may be rebinned to a 
lower resolution with only a small degradation of the results. 
For instance, rebinning to 10000 km/s (still centered on 642 and 
443 nm) gives p = 0.94 and cr(AyUfi) = 0.145 mag. Away from 
peak brightness the correlations weaken and the standardization 
slope y changes, but the 5(1 + z) day window used here is ad- 
equate for scheduling spectroscopy based on a pre-maximum 
lightcurve in any band. Although currently available high red- 
shift SN la spectra are insufficient to apply this method, fu- 
ture spectrophotometric instruments in space or on large ground- 
based telescopes could achieve these requirements and use this 
method to build an SN la Hubble diagram with the potential for 
significantly improved statistical power and systematics control. 



4. Conclusions 

We have performed a search for spectral flux ratios that corre- 
late with observed variability in SN la absolute magnitudes us- 
ing separate training and validation datasets totaling 58 homo- 
geneously observed SNe with 0.02 < z < 0.09. The results re- 
veal single parameter corrections that achieve results better than 
those from light curve shape and color parameters, even when 
only a single spectrum is used to measure and standardize the 
peak magnitude. The combination of a flux ratio and broad-band 
color provides even better results, significantly outperforming 
standard corrections for light curve shape and color. The best 
single flux ratio 7?642/443 results in a Hubble diagram residual 
scatter of 0.125 + 0.01 1 mag; the best combination of color plus 
a ratio provides 0.1 19 + 0.01 1 mag. 

These results highlight the power of spectrophotometric 
SN la observations. Compared to lightcurve shape and color cor- 
rections, our flux ratio standardization represents a factor of 1.8 
improvement in the statistical weight per SN, and a 35% im- 
provement in the upper bound on the potential bias from uncor- 
rected SN la variability. These results work well for both ex- 
trinsic dust extinction and intrinsic varations such as SN 1991T- 
like and SN 1999aa-like SNe, and they effectively correct SNe la 
over a range of light curve shapes and colors, host galaxy proper- 
ties, and intrinsic subtypes, offering the promise that this method 
will be similarly powerful for future high-z SN la observations. 
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Table 2. Supernovae used in this work, including their CMB-frame redshift z C mb, SALT2 fit parameters X\ and c, Hubble residuals A/i B , and flux ratios ft^/w 



Name 



ZCMB 



xi 



ft 



■642/443 



642/443 



Discovery 



Training SNe 

SN2005hc 

SN2006dm 

SN2007bd 

SN2007kk 

SNF20060908-004 

SNF2006 102 1-003 

SNF20061 11 1-002 

SNF20070424-003 

SNF20070531-011 

SNF20070630-006 

SNF200707 12-003 

SNF200707 17-003 

SNF20070727-016 

SNF20070803-005 
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Table 2. continued. 
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Notes: The Hubble residuals are before any corrections for c, x\, or "R. The peak magnitudes for the literature SNe Hubble diagram 
fits are from Hick en et al.l (120091) and lJha et al.l (120061) . which did not use SALT2. The ^42/443 Hubble diagram fits were performed 
separately for the literature and SNfactory observations; the reported A/jb residuals are with respect to the individual fits. The full 
table available from PDS includes all flux ratios from Table Q]for the SNe observed by the SNfactory. 



Table 2. continued. 



Name Zcmb x\ c A^s ^642/443 ^642/443 Discovery 

Discov erie s and redshifts are from SNfactory observations except as otherwise noted. Literature SNe redshifts are via Hicken et al. 
(120091) and lJha et al] (120061) . Other redshift citations via NED except for SDSS: 

"Theureau, G. et al. 1998, A&AS, 130, 333; 

fe Falco, E. E., et al. 1999, PASP, 111, 438; 

'Huchra, J. P., Vogeley, M. S., & Geller, M. J. 1999, ApJS, 121, 287; 
rf SDSS; 

e O'Neil, K. 2004, AJ, 128, 2080; 

^de Vaucouleurs, G. et al. 1992, VizieR Online Data Calalog, 7137, 0; 
SMathewson, D. S., & Ford, V. L. 1996, ApJS, 107, 97; 
''Keel, W. C. 1996, ApJS, 106, 27; 
'Jones, D. H. et al. 2005, PASA, 22, 277; 



